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Abstract. The slow and dense wind from a symbiotic red giant (RG) can be significantly deflected toward the orbital plane 
by the gravitational pull of the companion star In such an environment, the ionizing radiation from the companion creates a 
highly asymmetric Hii region. We present three-dimensional models of Hii regions in symbiotic S-type systems, for which 
we calculate radio maps and radio spectra. We show that the standard assumption of spherically symmetric RG wind results in 
wrong shapes, sizes and spectra of ionized regions, which in turn affects the observational estimates of orbital separation and 
mass loss rate. A sample of radio maps and radio spectra of our models is presented and the results are discussed in relation to 
observational data. 
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1. Introduction 

Symbiotic stars are binary systems consisting of a red giant 
(RG) as the primary, and a post-AGB star or MS-dwarf as the 
secondary. The RG primary is the source of a high density 
wind accreted by the secondary. In closer binaries (a < 10 AU) 
the wind is significantly deflected toward t he orbital plane by 
the gr avitational pull of the secondary (see Sawryszczak et al 
2002, and references therein). This process results in an en- 



hanced density region at the orbital plane, with a disk-like den- 
sity maximum surrounding the secondary. 

Radio emission from symbiotic stars is dominated by ra- 
diation from the ionized gas. The spectrum of an Hii region 
consists of three parts: optically thick at low frequencies (with 
spectral index a > 0.6), optically thin at high frequencies (with 
a ^ -0.1), and optically semi-thick at intermediate frequencies 
around at a turnover frequency Vt, with a smoothly decreasing 
from its optically thick value to -0.1. Taylor & Seaquist ( 1984) 



show that orbital separation a, ionizing luminosity Lph, and the 
ratio of mass-loss rate to wind velocity ^ may be derived from 
such observables as turnover frequency, monochromatic flux at 
that frequency ^t, and the value of the optically thick spectral 
index. 

However, their calculations are based on the spherical wind 



model originally developed by Seaquist et al. (1984) (hereafter 



STB), which, as we mentioned above, is incorrect in closer sys- 
tems. Good examples of such systems are symbiotic S-type 



giant and having 


orbital pe- 


Belczytiski et al. 


200C|, and 



references therein). 
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Recently, [Vlikolajewska & Ivison (2001 ) have obtained for 
the first time the spectral energy distribution (SED) at wave- 
lengths between 6 cm and 0.85 mm for one of the prototypical 
S-type systems, CI Cyg, during quiescence. Their data allowed 
to determine Vt, the optically thin ff emission measure, from 
which a lower limit to Lph and an estimate of a were calculated 
within the STB framework. Unfortunately, the comparison of 
these estimates with the known orbital and stellar parameters 
of CI Cyg poses a serious problem for the STB model. In par- 
ticular, a is overestimated by a factor of ~ 35, whereas Lph 
(whose value is known from optical/UV H i emission lines and 
He II /7/+j^continuum) is underestimated by a factor of ~ 20. 

The above inconsistency strongly indicates that the STB 
model needs revision. The most obvious improvements are re- 
lated to the aspheri city of the wind, as sugge sted both by nu- 
merical work (see [Gawryszczak et al 2002, and references 
therein) and observations. From the observational point of view 
there is little doubt about the origin of the ionized medium. 
The mm radio emission shows some correlation with the mid- 
IR flux, and the radio luminosity increases with the K - [12] 
colour. This indicates that warm dust is involved in the mass 
flow along with the ionized gas, and suggest that the red giant 
is the source of this material ( Mikolajewska et al. 2002a). Such 
a picture is also supported by the spectral analysis of symbiotic 
nebulae which showed CNO abundances similar to those ob- 



served in normal red giants (Nussbaumer et al, 1988) Optical 



and radio imaging of symbiotic stars reveal aspherical nebulae 
often with bipolar lobes and jet-like components ( porradi et al 



1999). It is hard to explain how such structures might develop 
if the RG wind were spherical. Additional, indirect evidence 
for the asphericity of the RG wind comes from emission line 
profiles. It is based on the stationary, blueshifted H„ absorp- 
tion originating close to the orbital plane and indicating that 
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the ionized region is probably bounded on all sides by a sig- 
nificant amount of dense neutral material (see the discussion in 



Gawryszczak et al, 2002). 

In the present paper we account for the asphericity of the 
RG as resulting from the interaction with the secondary. With 
the help of 3-D hydrodynamical models we demonstrate that 
Hii regions ionized by the secondary in S-type systems must 
indeed have shapes, sizes and spectra significantly different 
from those obtained by STB. A representative sample of radio 
spectra and radio maps of our models is also generated. 

In Sect. ^ we discuss the numerical techniques used in this 
work and details of the physical setup of our models. The mod- 
els are presented in Sect. ^ In Sect. ^ the results of our simula- 
tions are compared to observed symbiotic systems and a brief 
summary of our work is provided. 

2. Methods and assumptions 

The aim of the present paper is to study the properties of H ii 
regions in symbiotic S-type stars. In such systems the Hii re- 
gion develops in a dense RG wind whose density distribution 
is modified by the gravitational pull of the secondary. As a 
highly aspherical distribution may result even in cases when the 
wind is originally spherically symmetric, a multidimensional 
approach to the problem is necessary. 

We begin with 3-D hydrodynamical simulations of RG 
winds interacting with the companion to the mass-losing star 
The simulations are performed with the help of the Smoothed 
Particle Hydrodynamics (SPH) technique. An inertial coordi- 
nate system is used, with the origin at the mass center of the 
binary. 

Our implementation of the SPH is based on variable 
smoothing length with spline kernel introduced by Monaghan 



& Lattanzio (1985). The number of neighbouring particles is 



set to 40. For practical reasons (CPU time), the total num- 
ber of particles in the computational domain is limited to 10^. 
Whenever this limit is exceeded, those most distant from the 
system are removed from the domain. The distance beyond 
which removals occur (rjem) stabilizes within several orbital 
periods of the binary, and an almost stationary rotating density 
distribution is obt ained. The validity of such an approach was 
demonstrated by [Gawryszczak et aL (2002). With the above 
limit imposed on the number of particles, rrem does not exceed 
60 AU (see Table p. 

Within 0.1 AU from the secondary we remove all particles 
with velocities directed toward the secondary. This procedure 
is applied to avoid prohibitively short time-steps. 

There are good reasons to believe that RG winds are intrin- 
sically aspherical (see e.g. Frankowski & Tylenda 2001; Dorfi 
& Hofner 1996; [Reimers et al. 2000, and references therein). 
The asphericity may result from internal processes taking place 
in the star and/or presence of the secondary. However, since the 
theory is not yet advanced enough for precise quantitative pre- 
dictions, these eff'ects are not included in our models. In partic- 
ular, we assume that although the RG rotates synchronously it 
does not fill its Roche lobe, so that the mass flow through the 
inner Lagrangian point can be neglected. In effect, the wind is 
represented by a stream of particles launched at a constant rate 



from points distributed randomly and uniformly on the RG sur- 
face. The initial temperature of the wind is 3000 K, and the ini- 
tial velocity of wind particles is set to 1 km s"' ( note that, as the 
wind propagates across the grid its velocity increases; a mea- 
sure of the terminal wind velocity is introduced and explained 
in the following paragraphs). Both temperature and velocity of 
the wind are constant on the surface of the star. 

Further, we assume that the wind is composed of dust and 
ideal monoatomic gas which are dynamically coupled (i.e. they 
move at the same velocity). The details of wind acceleration 
process are not followed. Instead, we assume that the gravity 
of the red giant is balanced by radiation pressure on the dust, 
so that it does not appear explicitly in our equations. Thus, if 
the RG were single, the only force acting on the wind would be 
due to gas-pressure gradients, whereas in a binary the wind is 
additionally subject to the gravity of the companion. Such an 
approximation may seem rather crude; however Gawryszczak 
et al. (2002) demonstrated that it leads to terminal velocities 



comparable to those obtained by Winters et al. ( 2000 ) for winds 
in which radiation pressure nearly balances the gravity of the 
wind source. 

Finally, we neglect radiation pressure from the secondary ( 
note that in the case of a white dwarf,for any reasonable grain 
size the force due to radiation pressure on dust grains is much 
smaller than gravity). Self-gravity of the wind and any explicit 
radiative heating or cooling are also neglected. The effects of 
cooling are accounted for by a nearly isothermal equation of 
state employed for the wind (p = /cp'', with y = 1.01). A likely 
source of strong nonaxisymmetric effects is magnetic field, as 
the MHD collimation may be much more efficient than the 



gravitational focusing discussed here (see jarcia-Segura et al 



1999; Garcfa-Segura & Lopez 200C ; Garcfa-Segura et al. 2001 



and references therein). However, there is no information avail- 
able about magnetic fields in symbiotic systems and we feel 
that it is too early to incorporate them into models: the number 
of free parameters would simply grow too large. 

3. Results 

All binaries considered here consist of a 1.0 Mq red giant and 
a 0.6 Mq white dwarf. The orbital separation a varies between 
1 and 4AU, wh ich is representative of S-ty pe symbiotic sys- 
tems. Following Mikolajewska et al. (2002a) a mass loss rate of 
M - 10"^ Mq yr ' is adopted for the RG wind, while the ioniz- 
ing luminosity of the secondary Lph varies between 4 ■ lO'*'* and 
8 • 10"^^ phots"'. The latter range is compatible with estimates 



based on mm/submm observations (Ivison & Seaquist 1995 



Mikolajewska & Ivison 2001 ), but is a factor of ~10 smaller 
than estimates based on studies of optical/UV Hi and Hen 



emission lines and /7/+j^ continuum (Mikolajewska & Ivison 
2001| ; [Mikolajewska et all|2002a[ ). 

The choice of the lower estimate is dictated by the size of 
our hydrodynamical wind models (60 H- 100 AU in diameter), 
which at Lph S lO'** would be practically entirely ionized. The 
parameters of symbiotic binary models obtained in this work 
are listed in Table |l] (v^v in column 3 is a measure of the wind 
velocity; see Section for an explanation). The radius of the 
red giant is in all cases equal to 0.5 AU. A simple method of 
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scaling our results to combinations of M and Lph not covered 
by actual calculations is given in Section 

3.1. RG wind models 

For each binary in Table |I] we integrate the equations of hydro- 
dynamics until the density distribution of the RG wind relaxes 
to a stationary state. Density distributions of the relaxed mod- 
els are shown in (Fig. The general effect of the secondary is 
to deflect the RG wind toward the orbital plane, and as a result 
the relaxed density distributions are flattened. This flattening 
weakens with increasing orbital separation because of the in- 
creasing ratio of the wind velocity at the secondary's orbit to 
the orbital velocity of the secondary. However, all models de- 
velop a high density region in the form of a disk around the 
secondary. 

Table 1. List of models 



Model a[AU] 


Viiv [kms ' 


] nemLAU] 


Lph 


[phot s 


'] 




A 1.0 


25.0 


33 


7.3 


• 10"^ - 


7.3 


10*"' 


B 2.0 


17.0 


43 


8.0 


• 10*4 


8.0 




C 4.0 


17.0 


54 


4.0 


■ 10*^ ^ 


4.0 





All models also display a three-dimensional spiral struc- 
ture resulting from the shock wave excited by the motion of the 
secondary across the RG wind. The spirals are a very general 
feature of slow winds in binary systems, and they were also 
observed in simulations reported by [Vlastrodemos & Morris" 
( |1998| , |1999| ) or pawryszczak et al.| ( |2002| ). As one might ex- 
pect, in closer binaries the spiral is more tightly wound than in 
the wider ones. 

In the closest binary funnel-like density minima oriented 
perpendicularly to the orbit are also observed. Two effects con- 
tribute to the formation of these funnels. First, in a close binary 
the wind is emitted from a source moving with a high orbital 
velocity, and as such it carries appreciable angular momentum. 
Away from the orbital plane the centrifugal force associated 
with that momentum causes the wind gas to clear the vicinity 
of the polar axis. Secondly, and unlike that of the primary, the 
gravity of the secondary is not balanced by radiation pressure, 
and the wind flowing perpendicularly to the orbit is retarded. 

Below we argue that all these features (flattening, spirals 
and funnels) are important factors influencing the observational 
properties of the Hn regions created in the RG wind by the 
ionizing radiation from the secondary. 



3.2. Shapes and sizes of H ii regions 

Once the RG wind has been relaxed, its density distribution is 
smoothed by averaging over several tens of consecutive time 
steps. The aim of this procedure is to reduce the numerical 
noise resulting from the "grainy" nature of the SPH models. 
The smoothed density is mapped onto a spherical grid centered 




Fig. 2. 3D view of the H ii region boundary in model B with 
Lph = 2.4- 10"*^ Aregionof lOx lOx lOAUis shown. 

on the secondary, with grid points spaced uniformly in angu- 
lar coordinates and concentrated at the secondary in the radial 
coordinate. Then, the source of ionizing photons is turned on 
at the location of the secondary, and the boundary of the Hn 
region is found in a local Stromgren sphere approximation (i.e. 
Stromgren radii are calculated for all directions defined by grid 
points in angular coordinates). An example result of this pro- 
cedure is shown in Fig. ^. Separate treatment of ionization and 
dynamics is justified by the fact that the recombination time is 
one order of magnitude shorter than the orbital period at the 
edge of our computational domain, and more than 3 orders of 
magnitude near the secondary. 

In the approach developed by STB, the shape and size of 
the H II region are uniquely determined by the value of the pa- 
rameter X, which is a combination of the mass loss rate, wind 
velocity, orbital separation, and ionizing luminosity: 



X 



aL. 



ph 



(1) 



where arec and are recombination coefficient and wind ve- 
locity (which is constant in space and time), while n and wh 
have their standard meaning. 

While STB assume a wind of constant velocity, in our mod- 
els the flow accelerates with the distance from the RG, and the 
acceleration proceeds at a variable (location-dependent) rate. 
This is due partly to thermal pressure gradients, and partly to 
the centrifugal force associated with the angular momentum 
transferred to the wind from the secondary; see Gawryszczak 



et al. (2002). In order to calculate X from our models, we use 
a density-weighted velocity Vav, averaged over the whole solid 
angle (the averaging is performed near the edge of the compu- 
tational domain). For each binary a range of Hii regions with 
different X-values may then be easily generated by allowing 




Fig. 1. Density distributions of models A, B and C. The logarithmic density scale spans 5 orders of magnitude. In each frame a 
region of 40 x 40 AU is shown. The system rotates counterclockwise. 
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for variations in Lph (alternatively, M may be scaled by multi- 
plying the density of the RG wind in each grid cell by the same 
factor). 

To obtain the STB counterpart of our model, for the same 
M a density distribution of the spherical wind from the pri- 
mary is calculated, assuming - Vav An appropriate part 
of that distribution is fed into the spherical grid centered on 
the secondary, and the same Stromgren-routine is employed to 
find the boundary of H ii region. We checked that the boundary 
found this way agrees to within a few percent with the bound- 
ary derived from equation (1) of STB. 

The two sets of Hn regions are compared in Fig. |[ The 
upper part of each frame shows our model, and the lower one 
- the corresponding model obtained within the framework of 
STB (note that all displayed models are symmetric with respect 
to the orbital plane). 

A clear general trend may be observed: for wide binaries 
and intense ionizing fluxes (see lower right corner of Fig. ^ the 
shapes and sizes of our Hii regions converge to those of STB. 
On the other hand, for close binaries and/or low ionizing lumi- 
nosities, our H 11 regions become entirely different from those 
of STB: they are not axially symmetric, and their shapes may 
be fairly complicated (see Fig. ^ and the left column of Fig. |]). 
Another discrepancy concerns a large part of the space to the 
right of the secondary in Fig. ^ which in our models is screened 
from UV photons by the dense disk around the secondary. This 
neutral region, best visible in the upper right frames of Fig. ||, 
does not have its counterpart in STB models, and it only disap- 
pears when Lph is large enough to ionize the disk. 

Note also that our H ii regions are generally smaller than 
those of STB. This is due to the fact that the density of our RG 
wind is enhanced in the orbital plane compared to the spheri- 
cally symmetric wind with the same mass loss rate and veloc- 
ity. However, in the closest binary with the lowest Lph our H ii 
region becomes much larger than its STB counterpart. This is 
because the low-density funnels, which are very well developed 
in this case, are easily ionized, causing the H n region to evolve 
into a pair of semi-infinite jet-like lobes (upper left frame of 
Fig.|). 

3.3. Spectra 

The spherical grid centered on the secondary and containing 
the smoothed density distribution is also used to derive the ff 
spectra of our models. Monochromatic volume emissivity of 
the ionized gas i s calculated with the help of formulae taken 
from [Lang ( 198C ). Subsequently, a viewing direction is chosen 
and the emissivity is integrated along the line of sight until the 
optical depth of 5 is reached (or across the whole H ii region 
in optically thin cases). Observed monochromatic intensities 
are calculated for an assumed distance of 1 kpc between the 
observer and the radiation source. 

The resulting spectra of all models presented in Fig. |^ are 
displayed in Fig. ^ (note that in both figures the models are 
arranged in the same way, but the whole Fig. ^ is rotated coun- 
terclockwise by 90°). For comparison, the corresponding STB 
spectra are also shown. All spectra may easily be scaled to other 




Fig. 5. Scaling of the spectra. Compared are variants of model 
B with M = lO^^Moyr"', Lph = 2.4 ■ lO''^ (solid line), 
M = 3 ■ lO^^Moyr-i, Lph = 2.16-10**' (dashed line) and 
M = 3 • 10-^ Mq yv-\ Lph = 2.16 ■ lO** (dot-dashed fine). AH 
three H ii regions have the same shape and size. Their spectra 
are practically identical, except that they are shifted by approx- 
imately one order of magnitude in Sy and a factor of ~ 3 in 
V for every factor of ~ 3 in M accompanied by one order of 
magnitude increase in Lph. 



values of M and Lph; the corresponding procedure is illustrated 
in Fig |. 

As one might expect, in wider systems with high Lph our 
spectra almost entirely converge to those of STB, while they 
are significantly different in close systems with low Lph. In such 
systems our spectral index may be smaller or larger than that of 
STB; the same is true for turnover frequency and intensity. To 
make things even more complicated, some spectra exhibit more 
than one turnover. This is because in some models the distribu- 
tion of density along the line of sight cannot be approximated 
by a single power law. 

In general, the differences between the spectra of the same 
system viewed "pole-on" and along the line joining the stars 
(with the primary in front of the secondary) are larger in our 
models than in those of STB (slopes and turnover frequencies 
are however almost insensitive to such a change of viewing di- 
rection). On the other hand, our models predict weaker spectral 
variations with the orbital phase. This behaviour is explained 
by the fact that in our models the distinguished direction is per- 
pendicular to the orbital plane, whereas the H ii regions of STB 
are axially symmetric with respect to the line joining the stars. 

3.4. Radio maps 

Based on density distributions inside the ionized regions and 
the above described integration technique we derive maps 
of monochromatic ff emission from our models at 5 GHz. 
Depending on viewing angle and parameters of the binary, the 
maps display a rich variety of morphologies, of which a sam- 
ple is shown in Fig. ^. In some frames the outer boundary of 
the RG wind model is visible as a regular circular contour In 
such cases part of the emission is obviously not accounted for; 
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B 






C 





X = 0.2 



X = 0.6 



X = 2.0 



Fig. 3. A sample of H ii region shapes. H ii region boundary is plotted in a plane which is perpendicular to the orbit and contains 
the centers of both stars. In every frame our model (top) is compared to its STB counterpart (bottom) obtained for the same value 
of X. The columns illustrate the effect of increasing luminosity, (or, equivalently, of increasing X). In each frame a region of 
80 X 80 AU is shown; ionized regions are white. 



however it never amounts to more that a few per cent of the 
total emission at 5 GHz. Like in the case of spectra, the inten- 
sities indicated in the maps correspond to a distance of 1 kpc 
between the observer and the symbiotic binary. The following 
observations can be made: 

- Bipolar structure with clearly separated lobes is a standard 
for systems with X < I. 

- When the line of sight is inclined to the orbital plane, the 
lobes may differ both in shape and intensity. The differences 



are primarily due to the complicated shape of the effectively 
emitting (r ~ 1) layer in optically thick parts of the Hn 
region (an example of such region is shown in Fig. 
A significant degree of asymmetry may also be seen in 
pole-on systems (where only one lobe is visible). The pri- 
mary cause of this asymmetry is the 3-D spiral described in 



Sect. 3.2 



The bipolar structure shown by our maps could be easily re- 
solved by MERLIN after its upgrade is complete (e-MERLIN). 
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Fig. 6. Selected radiomaps at 5 GHz. Radio intensity is shown on a logarithmic scale, with solid contours plotted every order 
of magnitude (except for frames 5a-c, where a Unear scale is used). Every frame covers an area of 80 x 80 AU. Left column: 
side view (line of sight perpendicular to the line joining the stars in the orbital plane). Right column: pole-on view (line of sight 
perpendicular to the orbit). Middle column: intermediate case with line of sight inclined 45° to the orbital plane. The jagged 
appearance of some maps is an artifact due partly to limited resolution of the grid on which the boundary of the H n region is 
defined. Model parameters: la-c: a = 1, M = 10"^ Lph = 7.3 ■ lO'^, X = 0.2 2a-c: a = 1, M = lO"'', Lph = 2.2 ■ 10^^ X = 0.6 
3a-c: a = 2, M = 10"^ Lph = 8.0 • lO'*^ X = 2.0 4a-c: a = 4, M = 10"^ Lph = 4.0 • 10''^ X = 2.0 5a-c: a = 4, M = 
10"^ Lph = 4.0 ■ 10''^ X = 2.0. Average radio luminosity integrated over frames is 0.01, 0.10, 0.20, 0.35 and 0.05 mJy for rows 
1-5, respectively. 
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Even now, MERLIN has sufficient resolution (40 mas at 5 GHz 
where it is most sensitive) but its sensitivity is too low. 
However, e-MERLIN, with its new 15 and 22 GHz receivers 
will be able to perform imaging at a resolution of ~ 8 AU at 
1 kpc for sources as weak as a few of /iJy. 

4. Discussion 

We have demonstrated that the STB model is not applicable to 
systems with a ^ 4 and X ^ 2.ln such systems the H n region is 
not axially symmetric, and in most cases it is much smaller than 
its STB counterpart. This is because the basic STB assumption 
of a spherically symmetric RG wind is violated. Even if the 
wind is emitted uniformly from the whole surface of the red 
giant, it immediately loses its original symmetry due to grav- 
itational attraction of the secondary. As a result, a large-scale 
density enhancement develops in the orbital plane, with many 
smaller scale substructures (funnels, spirals and a compact cir- 
cumsecondary disk). Factors not accounted for in our simula- 
tions (RG wind originally enhanced in the orbital plane and/or 
outflow through the inner Lagrangian point) can only amplify 
this effect. Based on Fig. ^the following qualitative predictions 
can be made: 

- In most S-type systems large amounts of neutral material 
should be observed at the orbital plane. This is because the 
circumsecondary disk and the orbital density enhancement 
are much harder to ionize than the spherically symmetric 
RG wind. This prediction holds for the range of parameters 
investigated here (M ~ lO^^MQyr"'; Lph ^ 10'*'' phots"') 
as well as for M and Lph scaled as explained in Fig. |5} 

- Both the brightness and spectra of the radio continuum 
emission from edge-on systems should show little variabil- 
ity with the orbital phase. This is because the preferred di- 
rection (in which the ionizing quanta propagate most eas- 
ily) is perpendicular to the orbital plane. Note that in the 
STB model the preferred direction is the line joining the 
stars (or rather its extension beyond the secondary). 

- High resolution imaging of edge-on systems with lower Lph 
should reveal two maxima of radio emission. In some sys- 
tems traces of spiral structure may be detected. 

In fact, the first prediction has akeady been observation- 
ally verified. For many S-type systems there are strong indica- 
tions for low values of X (^ 1/3; ^iirset et al. 1991) and/or 
H II which ar e ioniz ation-bounded at least in the orbital plane 
( Ivison et al. ( 1994 ). In accordance with Fig. ^ such nebulae 



have a lso been observed for larger X's. Recently, [Quiroga et al 



(2002) have shown that in the case of AR Pav, the central ab- 



sorption in H„ is formed in a neutral portion of the RG wind 
which is strongly concentrated towards the orbital plane, (with 
the total hydrogen column density A^(H) S a few x 10^' cm"^). 

In systems with a ^ 2 and X ^ I our radio spectra signif- 
icantly diverge from those of STB. First, the location and ex- 
tent of the transition region between optically thin and optically 
thick parts of the spectrum are different. Second, the shape of 
the transition region is more complicated in our models. Third, 
our models radiate much less energy in the optically thick part 



of the spectrum than their STB counterparts (the difference can 
amount to two orders of magnitude). 

The latter result can be understood as a simple consequence 
of the smaller spatial extent of our H ii regions (the optically 
thick flux is roughly proportional to the projected surface area 
of the ionized part of the wind). In principle, the remaining 
predictions could be verified by spectroscopic observations. 
Unfortunately, observational data, especially at higher frequen- 
cies, are still too poor to discriminate between various the- 
oretical models. The broadest survey available, reported by 
[Vlikolajewska et ak (2002a ^) is based on mm/submm obser- 
vations of a sample of 20 S-type systems in quiescence. Radio 
emission from these systems is found to be optically thick at 
least up to ~ 1.3 mm (243 GHz = v, = turnover frequency). 
This is compatible both with the STB model, which for typical 
quiescent S-type systems predicts Vt ^ 300 GHz ( Seaquist et al 



1993), and with our simulations (see Fig. 



At present the only system which allows for a quantitative 
analysis of the spectrum is CI Cyg, where Vt ~ 30 GHz. For 
the STB model the outcome of such an analysis is rather un- 
favourable. In their framework the orbital separation may be 
estimated from 



a = 3001 



10"* K 



Vt 

GHz 



(mJy) ( 



kpc 



AU, 



(2) 



where Lg is the electron temperature of the ionized wind, and 
Si is the flux received at the turnover frequency. This estimate 
should be accurate within a factor of 2 for a range of X covering 
two orders of magnitude. However, in the case of CI Cyg it 
overestimates a by a factor of 35 (70 AU compared to ~ 2 AU 
obtained from the spectroscopic orbit solution. On the other 
hand, if one inserts a — 2 in ^ and calculates the turnover 
flux, then the derived 5t is 10^ times lower than the actually 
observed one. 

The STB framework also provides an estimate for Lph 



-^ph 



10^ V„ 



lO^K 

2 



GHz 



^ A, phots- 
mJy /\kpc ' 



(3) 



where 5v is the flux received in the optically thin part of the 
spectrum at frequency v [[vison & Seaquist (1995). However, 
for CI Cyg formula (^ yields a value 20 times lower than that 
obtained from optical/UV H i and He ii emission lines and bf+ff 
continuum. 

Other inconsistencies between the observed spectrum of CI 
Cyg and the STB model are discussed by Mikolajewska & 
Ivison (2001), who suggest that the STB approach might not 
apply to this system because it is likely that the red giant fills 
or nearly fills its Roche lobe, and the mass loss is concentrated 
in a stream flowing through the inner Lagrangian point. If this 
is really the case, then our models, which predict concentra- 
tion of the outflow in the equatorial plane and formation of the 
circumsecondary disk, should provide a better approximation 
of the observed spectrum than that obtained within the STB 
framework. 
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Unfortunately, in this particular case they do not perform 
much better A glance at Fig. Q indicates that in order to ob- 
tain Vt as low as ~ 30 GHz we would have to assume an un- 
acceptably low Lph and/or M much lower than lO^^MQyr"'. 
A likely solution to this problem is suggested by the spectrum 
for a - \ AU, X - 0.2 in Fig. ^ which is the only one show- 
ing an excess in the optically thick part. With increasing M and 
Lph the excess would shift to higher frequencies, mimicking 
the turnover observed in CI Cyg. The fact that the spectrum 
obtained for a - I AU offers a better fitting possibility for the 
system with a - 2 AU seems to indicate that our model of the 
latter underestimates the density of the circumsecondary disk 
and/or overestimates the density of the funnels. Such a con- 
clusion is compatible with arguments concerning the outflow 
through the inner Lagrangian point. 

The deciding tests of our results will become possible when 
high-resolution radio maps of S-type systems are obtained. Our 
results suggest that S-type systems will be interesting targets 
for e-MERLlN and ALMA. 

While the models presented here are closer to reality than 
those obtained within the STB approach, there is still a con- 
siderable room for improvements. We have only marginally re- 
solved the disk around the secondary, whose shape and extent 
may strongly influence shape and extent of the H ii region. To 
improve on that, future work should include effects of radiation 
transfer and viscosity, increase the mass-resolution of the SPH 
code and, accordingly, decrease its time-step. Note, however, 
that with radiation transfer and viscosity included, the simple 
scaling relation presented in Fig. |]would not work. This means 
that RG wind models extending to hundreds of AU would have 
to be obtained in order to generate H ii regions emitting as much 
energy in the optically thin regime as the observed ones. Such 
simulations would push the CPU requirements beyond cur- 
rently acceptable limits. When they are hopefully done in the 
future and a reliable thermal structure of the flow is obtained, 
the thermal ionization should be taken into account, which at 
least in some cases may significantly contribute to the radio 
spectra. 
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